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We describe an effective resonant interaction between two localized wave modes of different nature: 
a plasmon-polariton at a metal surface and a self-focusing beam (spatial soliton) in a non-linear 
dielectric medium. Propagating in the same direction, they represent an exotic coupled-waveguide 
system, where the resonant interaction is controlled by the soliton amplitude. This non-linear system 
manifests hybridized plasmon-soliton eigenmodes, mutual conversion, and non-adiabatic switching, 
which offer exciting opportunities for manipulation of plasmons via spatial solitons. 
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Introduction. — Plasmonics is an important, quickly 
developing area of modern physics which offers promising 
applications in nano-optics and electronics [1-4] . It deals 
with the so-called surface plasmon-polaritons [5], i.e. col- 
lective oscillations of the electromagnetic field and elec- 
trons which propagate along a metal-dielectric surface 
and decay exponentially away from the surface. 

Plasmons are characterized by frequency uj and prop- 
agation constant along the interface, kp — kp{uj) > k 
{k = ^/sqlu /c is the wave number in the dielectric medium 
with the permittivity Eq). Since kp > fc, plasmons can 
only interact resonantly with evanescent electromagnetic 
waves in the dielectric medium [1-4] . Accordingly, there 
are two main methods for excitations of plasmons: (i) 
via the evanescent wave generated at the total internal 
reflection [6] ; and (ii) via a periodic structure producing 
evanescent modes [4,7]. The first method can be modi- 
fied if there is a dielectric waveguiding layer parallel to 
the metal-dielectric interface. Propagating modes of the 
waveguide have evanescent tails outside the waveguide. 
They can interact resonantly with plasmons at the metal 
surface, that is used for diagnosis of layered dielectric 
films [8]. The plasmon channel can also be regarded as 
an effective waveguide, so that this configuration can be 
assimilated to the problem of two coupled waveguides. 

In this Letter, we propose a new way of resonant in- 
teraction of plasmons with electromagnetic waves. In- 
stead of exploiting various inhomogeneities to generate 
evanescent modes, we introduce a non-linear dielectric 
medium. A nonlinear dielectric admits self-focusing so- 
lutions (spatial solitons) [9] which may propagate paral- 
lel to the metal-dielectric interface and, similarly to the 
modes of the dielectric waveguide, also have evanescent 
tails and propagation constant kg > k [10]. A 'quasi- 
waveguide' is formed by the soliton profile which modifies 
the effective dielectric constant of the non-linear medium. 
Hence, the soliton-plasmon configuration couples effec- 
tive linear and non-linear waveguides, as shown in Fig. 
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FIG. 1: (Color online.) Plasmon-soliton system as two cou- 
pled waveguides. 



1, and an interaction between the plasmon and soliton 
may occur at certain resonant parameters. In contrast to 
previous studies in non-linear plasmonics, where effects 
of nonlinearity on plasmon modes have been considered 
(see, e.g., [11]), here we examine interaction between two 
spatially-separated modes: soliton and linear plasmon. 
This system is essentially self-influencing - the coupling 
is controlled by the soliton amplitude rather than by the 
conflguration parameters. 

Model. — To create a model describing the plasmon- 
soliton system, Fig. 1, we adopt several simplifying as- 
sumptions. First, the waves propagate along the z coor- 
dinate and the wave electric fleld lies in the {x, z) plane, 
so that the y coordinate can be eliminated from further 
consideration. Second, we assume that the dielectric non- 
linearity is localized around a certain distance x — d from 
the metal surface a: = 0, so that it practically does not 
affect the plasmon field. The distance d is assumed to be 
larger than the characteristic widths of the plasmon and 
soliton fields, which guarantees an exponentially-small 
overlapping of their tails. Hence, the coupling is weak 
and can be treated perturbatively. 

Uncoupled plasmon and soliton fields, '^p and \E's 
[12], can be represented as '^p{x,z) = Cp{z)'ipp{x) 
and 5's(x,z) = Cs{z)'4!s(x,\cs\), where Cp^s are the z- 
dependent amplitudes and 'i/'p,s are the transverse pro- 
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files of the fields with the normalizations "0^(0) = 1 
and i^sid) = 1. The transverse profile of the plas- 
mon, ipp, represents two exponents decaying away from 
the metal-dielectric interface [1,2] {ipp — exp(— Kpx), 

Kp = ^Jk"^ — k'^ in the dielectric, at a; > 0), whereas 

the soliton profile is given by [10] ips — sech [ks{x — d)], 
Kg = k^y j/2\cs\, Fig. 1. Here 7 is the parameter of 
nonlinearity of the medium [13]. We emphasize that the 
transverse profile of the soliton depends on its amplitude 
|cs|. In what follows we assume that the amplitude varies 
slowly enough, so that one can use quasi-stationary adia- 
batic approximation considering \cs\ as a local parameter 
of the problem. 

Thus, the total wave field in the problem, ^P, can be 
represented via the ansatz 



^{X,z) = Cpiz)'lljp{x) + Cs{z)'lpsix, \Cs\) 



(1) 



In the zero approximation, when the plasmon-soliton 
coupling is negligible, amplitudes Cp and Cg are indepen- 
dent and obey, respectively, linear and non-linear oscil- 
lator equations. In the first order approximation, the 
plasmon and soliton fields become linearly coupled due 
to the spatial overlapping of '^p and and the trans- 
verse inhomogeneity of the medium along x (i.e., the con- 
trast of dielectric indices at the metal-dielectric interface, 
£m 7^ £q)- This results in the coupled oscillator equations 
for the amplitudes Cp and Cg: 

c'; + (3lcp^q{\cs\)cs , c': + f3U\cs\)cs = q{\cs\)cp . (2) 

Here the prime stands for the derivative with respect to 
the dimensionless coordinate — kz, Pp — kp/k > 1 and 
(3s = ks/k ~ 1 -I- 7|cs|^/4 (we assume that 7|cs|^ ^ 1) 
are the plasmon and soliton dimensionless propagation 
constants, and q <C 1 is the coupling coefficient. 

To determine the coupling coefficient q, note that 
the right-hand side of the first Eq. (2), i.e. gdc^Dcs, 
represents an external source that excites plasmons at 
the metal surface. Hence, up to a numerical factor, it 
should be equal to the soliton field at the metal surface, 
vj/^l^j^g = Csips\x=o- This yields the estimation 



q{\cs\) ^ iIJs\x=o ^ exp I 



(3) 



which is adopted below. Note that the weak-coupling 
approximation fails at small soliton amplitudes. 

Equations (l)-(3) represent the model describing the 
plasmon-soliton interaction. These equations are remi- 
niscent of the system of two weakly coupled nonlinear 
waveguides [14]. However, there are two essential pecu- 
liarities in the plasmon-soliton equations: (i) only one 
subsystem is nonlinear; (ii) the coupling coefficient de- 
pends on the soliton amplitude. 

An effective plasmon-soliton interaction and energy ex- 
change occurs only near the resonance /3p = Ps which is 
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FIG. 2: (Color online.) Collective modes in the plasmon- 
soliton system vs. the dimensionless soliton amplitude, (a) - 
Eigenvalues Eq. (7) (bold curves), eigenvalues of uncou- 
pled modes and Va (dashed lines), (b) - Amplitude ratios 
a± Eq. (8). (c) and (d) - Profiles of the total field |*(2;)|^ 
Eq. (1), corresponding to the points A-F at curves in (a) 
and (b). The plasmon tails at a; < are not depicted, since 
they are determined by the properties of a particular metal. 
Parameters are z/p = 0.2 and fed = 6 (^ros — 0.02). 



achieved at the soliton amplitude | Cs [res = •\/4(/3p — l)/7. 
Solitons with significantly larger or smaller amplitude are 
uncoupled from the plasmon. In the vicinity of resonance, 
/3s| ^ Eqs. (2) can be simpfified. Making substi- 
tution Cp^s(C) = Cp.s(C) exp(iC) and assuming that new 



amplitudes Cp,s vary slowly (|Cp ,,[ <C |Cp_, 
be neglected), we arrive at the following equations: 



Cp 



Vp -q{\Cs\)l2 
-q{\Cs\)/2 Vs{\Cs\) 



Cp 

Cs 



(4) 



Here 



/3p - K 1 



-l=7|a|V4<l (5) 



are the small deviations of the dimensionless propagation 
constants from unity. The first inequality (5) implies 
plasmons in the long- wave region of their spectrum close 
to the light cone, whereas the second inequality indicates 
the weakness of nonlinearity . 

Equation (4) has the form of a vector non-diffractive 
non-linear Schrodinger equation with a non-diagonal 
Hamiltonian typical in quantum two- level systems [15]. 
It possesses the integral of motion |Cp|^ -I- |Cs|^ = const, 
which is associated with the conservation of the total en- 
ergy under the plasmon-soliton interaction. 

Eigenmodes. — Modes of Eq. (4) are obtained via 
substitution Cp^siQ) = ^p,sexp(wC). This yields the 
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characteristic equation determining the eigenvalues, 



--i\Cs\) 



2 

(6) 

and the ratios of the plasmon and soliton amplitudes de- 
termining the eigenvectors, 



--i\c.\)^^ = - 



q{\Cs\) 



2[z/±(|C.|)-z/p 



(7) 



Unlike linear systems, here the eigenvectors cannot be 
normalized arbitrarily because the right-hand side of Eq. 
(7) depends on the amplitude \Cs\ = \As\. 

Equations (6) and (7) describe the collective (hy- 
bridized) plasmon-soliton modes appearing due to the 
coupling. They have a standard form typical for a linear 
two-level problem with close eigenvalues and small cou- 
pling [15]. In our case, the problem is essentially non- 
linear and the properties of the solution depend on the 
soliton amplitude \Cs\ — |^s|- In other words, the soliton 
amplitude plays the role of the driving parameter. 

Figure 2 shows the eigenvalues v^, Eq. (6), together 
with Vp and z/g, amplitude ratios a^, Eq. (7), and field 
profiles |\l/(a;)p, Eq. (1), of the collective plasmon-soliton 
modes as dependent on the soliton amplitude \Cs\. In the 
interaction region, |fp — J^s| ^ Qres, the eigenvalues avoid 
crossing (the minimal frequency gap grcs = 9(1 Colics) = 



exp (- 



is achieved at resonance), and eigen- 



modes are collective excitations with \Ap\ ~ At 
resonance, I's — Vp, the two modes show their closest 
intensity distributions (B and E) corresponding to the 
symmetric and antisymmetric combination of plasmon 
and soliton states with equal amplitudes. Away from 



the interaction region. 



^ (^res, the eigenvalues 



V tend to the asymptotes Vp and Vs representing pure 
plasmon and soliton states. Accordingly, the energy is 
concentrated basically in either the plasmon (D, C) or 
the soliton (A, F) channel. When passing through the 
interaction region along the same branch in Fig. la, the 
eigemodes exchange the amplitudes. Fig. 2b, and a near- 
soliton mode transmutes into the plasmon one and vice 
versa. This offers various possibilities for energy conver- 
sion between the plasmon and soliton channels. 

Dynamics and dissipation. — First, we consider the 
case when the soliton channel is excited at the input with 
a near-resonance amplitude Cs(0): 



Cp(0) = 0, |C,(0)- ICIresI « |C,|, 



(8) 



Pure soliton is not an eigenmode, and initial conditions 
(8) excite a mixture of and "— " eigenmodes with 
close frequencies. In a linear problem this would lead 
to a superposition of two modes and a harmonic low- 
frequency beating. But in the plasmon-soliton system, 
the superposition principle is invalid, and the solution of 
Eq. (4) with Eq. (8) represent a self-modulated non- 
linear beating causing a periodic mutual conversion of 




FIG. 3: (Color online.) Nonlinear beating of the plasmon 
(lower curves) and soliton (upper curves) dimensionless am- 
plitudes after the excitation of the soliton channel near reso- 
nance. Eqs. (4) and (8) are solved numerically with Vp — 0.2 
and values of 7|Cs(0)p/4 as indicated in the boxes. Values of 
kd correspond to grcs — 0.02 (a) and grcs — 0.003 (b). 



energy between the soliton and plasmon channels, Fig. 3. 
The maximal conversion amplitude is reached at a certain 
value Cs (O)niax > |Cs|res- K corresponds to the maximal 
beating period ~ 27r/gres and it drops abruptly for 
higher Cs(0). The beating amplitude decreases with the 
resonant coupling ^i-cs decrease. Fig. 3a and b. 

Second, soliton-plasmon conversion occurs under ma- 
nipulations with the soliton amplitude - the driving pa- 
rameter of the system. One could expect that slow 
changes of the soliton amplitude will result in an adia- 
batic transformations of the eigenmodes along the "+" or 
"— " dispersion curves. Fig. 2a, and the soliton will meta- 
morphose into plasmon when passing over the resonance. 
For instance, the soliton is excited with an off-resonance 
amplitude corresponding to the zone A in Fig. 2, i.e. 



Cp(0)=0, CsiO) ~ \CsU ^ \Cs\r 



(9) 



These initial conditions approximately match the "+" 
eigenmode. Then, we introduce a small absorption in 
the medium which makes the soliton amplitude slowly 
decreasing along \Cs\ — \Cs\{C)- The losses in the plas- 
mon and soliton channels are described by adding small 
imaginary parts fjp.s to their propagation constants: 



^p.s 



-'p.s J 



(10) 



The weak dissipation (10) does not affect the local char- 
acteristics of the waves, which change adiabatically with 
|Cs|(C) according to the eigenmodes (6) and (7). 

With Eq. (10), this leads to the transformation of 
states as A ^ B ^ C, Fig. 2, i.e. almost 100% 
conversion of energy from the soliton to the plasmon 
channel. This would be a perfect switch, but the non- 
adiabatic Landau-Zener transitions between the and 
"— " eigenmodes may appear in the vicinity of resonance, 
|fp — t^sl ^ (/res [15,16]. The probability of the Landau- 
Zener transitions can be estimated using the parameter 
of adiabaticity /i [17]: 



,±'i 



(11) 
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FIG. 4: (Color online.) Adiabatic dynamics in the plasmon- 
soliton system with small dissipation, Eqs. (4), (9)-(ll), — 
0.2 and kd = 6. Initial soliton (9) attenuates slowly until it 
reaches the resonant amplitude. Then, only a part of the 
energy is abruptly transferred to the plasmon channel (/ires ^ 
1) with a subsequent beating in the mix of "+" and "— " 
eigenmodes, cf. Fig. 3. The value of the adiabatic parameter 
(thin dotted line) indicates the effectiveness of the energy 
exchange between the eigenmodes. 



chalcogenide glass, eo 5 ~ 10). In this range, plas- 
mons are characterized by the propagation constant close 
to unity, ^ 10^'^ and relatively small dissipation 
(Tp ~ 0.1. Solitons with i^s = 7|csp/4 ^ 10"'^ can be 
produced with typical Kerr nonlinearities. Transverse 
evanescent tails of plasmon and soliton are determined 
by Kp ^ Ks ^ 0.1/im~^, and propagation of the soli- 
ton at the distance d ~ 20 40/xm {kd ~ 50 100) 
from the metal surface yields the coupling coefficient 
Qrcs ~ 10~^ 10~^. Although the plasmon dissipation is 
large enough, ap ^ gres [18], we have found, by numerical 
simulations of Eqs. (4) with (10), that our conclusions on 
the plasmon-soliton coupling and dynamics (e.g.. Fig. 4b) 
remain qualitatively true at these realistic parameters. 

This work was supported by contracts FIS2005-01189 
and TIN2006-12890 from the Government of Spain, and 
by the Australian Research Council. 



If <C 1, the "+" and "— " eigenmodes evolve indepen- 
dently according to the adiabatic evolution. On the con- 
trary, if /i 3> 1, the wave undergoes diabatic evolution 
and practically all the energy stored in the eigen- 
mode is converted into the "— " one when passing the 
resonance. This means the transformation of states as 
A B + E ^ F, Fig. 2, i.e. the non-adiabatic tran- 
sition keeps the energy in the soliton channel. In the 
intermediate case, /i ^ 1, the output wave will present a 
mixture of the "+" and "— " modes: A ~* C + F. 

Surprisingly, the adiabatic parameter (11) cannot be 
made small in the plasmon-soliton system, no matter how 
weak the dissipation is. The point is that the soliton am- 
plitude I Cs I dramatically changes in the vicinity of reso- 
nance due to soliton-to-plasmon transmutation along the 
eigenmode, and the adiabatic regime is never achieved. 
This is a purely nonlinear effect - the driving parameter 
is a part of the solution. Numerical simulations of Eq. (4) 
with Eqs. (9)-(ll) confirm this conclusion. Fig. 4. In the 
vicinity of resonance, /iics ^ 1 which results in a par- 
tial soliton-to-plasmon switching and the Landau-Zener 
transition with the consequent nonlinear beating between 
and "— " eigemodes. Smaller values of the coupling 
gres result in higher /ij-es and smaller part of energy trans- 
ferred to the plasmon. 

Discussion. — To summarize, we have shown the 
possibility of resonant interaction between a plasmon- 
polariton at a metal surface and a parallel self-focusing 
beam in a nonlinear dielectric. A simple two-level model 
reveals hybridized plasmon-soliton eigenmodes and their 
complex non-linear dynamics which offers plasmon exci- 
tation and control using spatial solitons. 

An effective soliton-to-plasmon coupling and energy 
conversion can be achieved in the far IR frequency range 
(the wavelength in vacuum is Aq ~ 5 lO^m), by us- 
ing a gold interface with a non- linear dielectric (e.g., a 
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